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Abstract 



We present results of numerical simulations to estimate scaling exponents 
associated with driven surface growth in two spatial dimensions. We have 
simulated the restricted solid-on-solid growth model and used the time and 
system size dependent interface width, and the equal time height correlation 
function to determine the exponents. We also discuss the influence of vari- 
ous functional fitting ansatzes to the correlation function. Our best estimates 
agree with the results of Forrest and Tang obtained for a different growth 
model. 



PACS numbers: 05.40, 05.70L, 61.50C, 64.60H. 
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Kinetic roughening of interfaces is an ubiquitos phenomenon which takes 
place under a variety of non - equihbrium conditions [0. One of the most 
commonly used examples of a class of systems undergoing this process is 
driven surface growth far from equilibrium, where the surface current is non- 



In this case the relevant mapping of the problem in the continuum limit is 
given by the Kardar-Parisi-Zhang (KPZ) equation 



where v and A are constants, and r/ is a random noise term with (//(r, t)ri[f', t')) = 
2D6'^{f— f ')6{t — t'). The height variable h{f,t) describing the surface is a 
function of time and two-dimensional vector r, and the total dimensionality 
is denoted hj d = 2 + 1. 

Due to the scaling relation for the interface width 0] 



where L is the system size, and the scaling relation z+x = 2, there is only one 
independent scaling exponent for the problem. It is usually determined either 
from w{t) ~ t^, where (3 = x/ z or from the steady-state limit w{L) ~ L^. 

The scaling exponents are exactly known only in d = 1 + 1 dimensions, where 
/3 = 1/3 0. A lot of analytic and numerical work has been carried out in 
order to establish general, dimension dependent values for the exponents |[l[]. 
Monte Carlo simulations of discrete growth models have proven useful for 
this purpose. For example, extensive work [^] on the hypercubic stacking 



conserved and diffusion rates are very slow compared to the driving force . 




(1) 



w{L,t) ~ L"" 



f{tL-^), 



(2) 
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model has given (3{3) = 0.240(1), and (3{A) = 0.180(5). Ala-Nissila et al. [| 
simulated the restricted solid-on-solid growth (GRSOS) model up to d = 7+ 
1. By concentrating on > 3 + 1, they obtained /9(4) = 0.180(2) in excellent 
agreement with Ref. 0, but not with the conjecture (3{d) = l/{d+l) 
and showed that there is no upper critical dimension up to (i = 7 + 1. Their 
high accuracy data for (i = 3 + 1 was based on a novel fitting ansatz for the 
equal time correlation function 

G{r,t)^{[h{r' + r,t)-h{r',t)f)r>, (3) 

averaged over r', as 

Gi{r,t) = al(^){tanh[6l(^)r^lW]}^^ (4) 

where ai(t), bi{t) and 7i(t) = 2xi(t)/xi are fitting parameters, and Xi is 
fixed. This functional form gives, after fixing xi, an estimate of /5 through 
ai(t) ~ t^^, and also for x ^ Xi ^ 0- 

In this note, our purpose is to present new simulation data for the GRSOS 
model in two spatial dimensions, and estimate the corresponding scaling ex- 
ponents. This case is particularly interesting for its potential for experimental 
realizations, and also from a theoretical point of view. In contrast to other 
recent work claiming j3 = 0.25 our best estimate /5(3) = 0.240(2) is in 
excellent agreement with the hypercubic stacking model |^ , although finite- 
size effects seem to be somewhat pronounced. We also discuss the influence 
of different forms of fitting functions, in addition to Eq. (4), to exponents 
extracted from the time-dependent correlation function. 
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First, we used the time dependent width w'^{t) for several system sizes to 
determine (3. The results of least -squares fitting are presented in Table 1. 
In Fig. 1 we show these values plotted against 1/L. Result for the largest 
system size studied L = 2000 comes already very close to 0.240 as obtained in 
Ref . P , although statistical errors increase considerably for largest systems. 
We note that analysis of the data in the form of log[w^(2t) —w'^{t)] vs. log(t) 
as in Ref. failed to produce any consistent results. 

To corroborate these findings, we next calculated the saturated width w'^{L) ~ 
for L = 125, 250, and 500. These data give u;^ = 5.15(2), 8.67(2), and 
15.1(8), where the errors have been estimated from fluctuations between con- 
secutive runs. From a least squares fit we obtain x(3) = 0.387(2), which gives 
/3(3) = 0.240(2), in complete accordance with the time dependent width. 

As previously shown , the fitting ansatz (4) can be used to obtain accurate 
estimates of f3 even for relatively small systems. In Ref. 0, accurate results 
for f3 were obtained by using Eq. (4). It was also shown that the values of Xi 
obtained were somewhat smaller than those corresponding to /3 (as extracted 
from ai{t)), except in = 1 + 1 dimensions. In the previous work, xi was 
fixed to be Xi = 1 (c? = 1 + 1) or 1/2 (c? > 3 + 1). In this work, we let Xi 
vary and also extend the original fitting ansatz to include the following new 
fitting functions: 

G2{r,t) = a2(t){l -exp[-62(t)r^^(*^]r, (5) 

and 
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G3{r,t) = as{t){-- + arctan[exp(-63(^)r^^(*))]}"^ (6) 

where a2(t),a3{t),b2{t), 63 (t), 72 (t) = 2x2{t)/x2, and 73 (t) = 2x3it)/x3 are 
new fitting parameters. To perform tlie fitting, we calculated 3000 averages 
of the correlation function (3) for L = 100, and 540 averages for L = 200 
and 500. 

To test the quality of the fitting functions, we first fixed Xi = 0.387 for 
each function, and allowed Xj's to vary. Data for L = 200 was used. For 
each function then, was fixed corresponding to an average value obtained 
between 248 and 600 Monte Carlo time steps per site (MCS/s), where typical 
variations of less than about 10 %. This gives xi = 0.7389, X2 = 

0.4781, and X3 = 0.6484, which values were in turn used to obtain xi = 
0.387(8), X2 = 0.387(5), and xs = 0.386(8), correspondingly, as average 
values over 248-800 MCS/s. Next, we calculated estimates for the exponent 
z averaging results for c = 0.9 and 0.95 over 248-800 MCS/s (as explained in 
Ref. II), and obtained z = 1.60(2), 1.62(2), and 1.61(2) for fitting functions 

(4) , (5), and (6), respectively. These values obey the scaling relation x + z, 
giving 1.99(3), 2.01(3), and 2.00(3), respectively. This demonstrates that the 
quality of each fitting function is very good. 

Next, we fixed X2 as above and calculated 02 (t) for the fitting funtion of Eq. 

(5) , which gave the smallest variation for x. Results for the other functions 
are consistent. For L = 100 and 200, both the original correlation function 
and the fitting data are very smooth, and least squares fitting to 02 (t) gives 
/3 = 0.242(2) and 0.240(2), respectively. For L = 500, the data are not as 
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good, but fitting to a relatively straight region of the log-log curve gives 
0.238(1) (the error bars are purely statistical). These values are in excellent 
agreement with data obtained from the width for the largest systems in Table 
1. 

As a final check, we calculated the scaling function of Eq. (2) as shown in 
Figs. 2(a) and (b). For the range of system sizes studied in Table 1, (3 = 0.24 
in Fig. 2(a) gave clearly better scaling than f3 = 0.25 of Fig. 2(b). In the 
inset of each figure, we also show the scaling functions for L = 1000 and 2000 
only. For these two largest system sizes, the accuracy of the data does not 
allow us to distinguish between the two values of the exponent. 

To summarize, we have presented results of rather extensive simulations of 
the GRSOS model for system sizes up to L = 2000. Our best estimate 
/9(3) = 0.240(2) comes out from various independent ways of determining the 
scaling exponents, and is in excellent agreement with Ref. |^ for the range of 
system sizes considered in the present work. Unfortunately, a straightforward 
extension of this work to larger systems becomes prohibitive, since several 
hundreds of hours of computer time in RISC workstations was required here 
already. 

Acknowledgements: This work has been partly supported by the Academy 
of Finland. Centre for Scientific Computing Co. and Tampere University of 
Technology are acknowledged for the computer resources. We wish to thank 
the referee for useful suggestions, and the Department of Physics at McGill 
University for hospitality while this work was being completed. 

Note added in proof: A single run for a. L = 4000 system gives results fully 
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consistent with L = 2000, with (3 increasing at later times accompanied by 
very large fluctuations in w{t). 
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Figure Captions 



Fig. 1. Estimates of /3(3) as obtained from least-squares fits to w'^{t) (see 
Table 1 for details). 

Fig. 2(a). Scaling of w(L,t) for L = 100, 200, 300, 1000, and 2000, with 
P = 0.24, and (b) the same data for /3 = 0.25. Insets show scaling between 
L = 1000 and 2000 only. 
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Table Caption 



Table 1. Results of least-squares fitting to the time and system-size depen- 
dent width w'^{L,t). Error bars are purely statistical. 



8 



References 



[1] J. Krug and H. Spohn, in Solids Far From Equilibrium: Growth, Mor- 
phology and Defects, C. Godreche ed. (Cambridge University Press, 
Cambridge 1991); P. Meakin, Phys. Reps. 235, 191 (1993). 

[2] J. Villain, J. Phys. (Paris) II, 19 (1991). 

[3] D. Forster, D. R. Nelson, and J. M. Stephen, Phys. Rev. A 16, 732 
(1977); M. Kardar, G. Parisi, and Y. C. Zhang, Phys. Rev. Lett. 56, 889 
(1986); E. Medina, T. Hwa, M. Kardar, and Y. C. Zhang, Phys. Rev. A 
39, 3053 (1989). 

[4] F. Family and T. Vicsek, J. Phys. A 18, L57 (1985). 

[5] B. M. Forrest and L-H. Tang, Phys. Rev. Lett. 64, 1405 (1990); L-H. 
Tang, B. M. Forrest, and D. E. Wolf, Phys. Rev. A 45, 7162 (1992). 

[6] T. Ala-Nissila, T. Hjelt, and J. M. Kosterlitz, Europhys. Lett. 19 (1), 1 
(1993); T. Ala-Nissila, T. Hjelt, J. M. Kosterlitz, and O. Venalainen, J. 
Stat. Phys. 72, 207 (1993). 

[7] J. M. Kim and J. M. Kosterlitz, Phys. Rev. Lett. 62, 2289 (1989); J. M. 
Kim, J. M. Kosterlitz, and T. Ala-Nissila, J. Phys. A 24, 5569 (1991). 

[8] J. M. Kim, M. A. Moore, and A. J. Bray, Phys. Rev. A 44, 2345 (1991); 
D. Ko and F. Seno, cond-matObabbage . sissa. it no. 9312054 (1993). 



9 



L 


m 


Number of runs 


100 


0.226(1) 


3000 


200 


0.231(2) 


2750 


300 


0.232(2) 


400 


1000 


0.236(2) 


25 


2000 


0.239(3) 
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Table 1: 
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